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QO ' Abstract 

We present a new parton model approach for nuclear collisions at RHIC energies (and beyond). It 
is a selfconsistent treatment using the same formalism for calculating cross sections like (Jtot or (Jinci 
►^ ' and, on the other hand, particle production. Actually, the latter one is based on an expression for the 

QQ I total cross section, expanded in terms of cut Feynman diagrams. Dominant diagrams are assumed to 

^SJ . be composed of parton ladders between any pair of nucleons, with ordered virtualities from both ends 

^^ ' of the ladder. 

^' 

o 

ON ■ 1 Introduction 

C~| ' The standard parton model approach to proton-proton or also nucleaus-nucleus scattering amounts to 

CL|| presenting the partons of projectile and target by momentum distribution functions, /^ and /s, and 

f^ • calculating inclusive cross sections as |sjo82, dur87|| 
(D ■ 

^ ■ rr:.... ^\^ I dt I d^+ I d'r-f\(T+ 0'^)fl(T~ 0^\^!f^(T+ 
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rS I where da-ij/dt is the elementary parton-parton cross section. This simple factorisation formula is the result 

j^ ■ of cancellations of compHcated diagrams (AGK cancellations) and hides therefore the complicated multiple 

scattering structure of the reaction. The most obvious manifestion of such a structure is the fact that the 
inclusive cross section exceeds the total one, so the average number of elementary interactions must be big- 
ger than one. The usual solution is the so-called eikonaHzation, which amounts to re-introducing multiple 
scattering, however, based on the above formula for the inclusive cross section. The latter formula has the 
disdvantage of having lost many important details about the multiple scattering aspects. Constructing in 
this way a model (in particular an event generator) for particle production is completely arbitrary. For 
example, energy obviously must be conserved, but the QCD formulas will not provide the slightest hint 
how to realize energy conservation in a partic ular m u ltiple s cattering event. 

This problem has first been discussed in | abr92| ,[ bra9C |. The authors claim that following from the 



nonplanar structure of the corresponding diagrams, conserving energy and momentum in a consistent way 
is crucial, and therefore the incident energy has to be shared between the different elementary interactions, 
both real and virtual ones. 

Following these ideas, we provide in this paper a rigorous treatment of the multiple scattering aspect, 
such that questions as energy conservation are clearly determined by the rules of field theory, removing 
the arbitrariness of the procedures applied so far. The general idea is as follows: starting point is an 
expression for the inelastic cross section for nucleus-nucleus scattering at very high energies, expressed in 
terms of cut Feynman diagrams. Clearly at this point certain assumptions concerning the dominace of 



certain classes of diagrams is employed. But this assumptions really define the model, what follows is just 
the application of the rules of field theory, calculating diagrams, making partial summations, and finally 
interpreting these partial sums as "partial cross sections" fo r certa in physical processes. Production of 
particles (on the parton level) is thus completely determined | wer97| . 



2 The Elementary Cut Diagram: the Parton Ladder 

We want to write down an expression for the inelasic cross secion in nucleus-nucleus (including nucleon- 
nucleon) scattering in terms of cut Feynman diagrams. We will assume that dominant contributions will 
come from certain classes of diagrams, which are composed of so-called "elementary diagrams", also referred 
to as "elementary interactions". We will therefore first dicuss the elementary diagrams, before introducing 
the multiple scattering theory. A complete diagram will certainly contain cut and uncut elementary 
diagrams. We will investigate explicitely only cut diagrams and use a relation between cut and uncut 
diagrams to calculate the latter ones. 

We assume an elementary interaction to be represented by a parton ladder with "soft ends", see fig 
|l|. The central part is a parton ladder with ordered virtualities, such that the highest virtuality is at 
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Figure 1: The elementary cut diagram: parton ladder plus "soft ends". 



the center and the virtualities are decreasing towards the ends of the ladder. This part of the diagram 
can be calculated using perturbative techniques of QCD. Since the virtualities are decreasing towards the 
ends, one reaches finally values where perturbative calculation can no longer be employed, although the 
longitudinal momentum fraction of the corresponding parton may be much smaller than one. This means 
there is still a large mass "object" between the first parton of the ladder and the nucleon, however, with 
small virtualities involved |lan94|. The most naturel candidate for such an object is the soft Pomeron, 
which can not be calculated from first principles, but where reasonable parametrisations exist, based on 
general considerations of scattering matrices in the limit of very high energies. So, the mathematical 
expression corresponding to the cut diagram of fig. |^ is 



i?.s, 



{s,x\x-,b) = fdb'j2f^ f^Eu ( ^,b')EU ( ^,&-&>;iddc.(4^r-), (1 



where E'soft represents the soft Pomeron at each end of the ladder and CTj'j^ddcr the parton ladder itself, the 
precise definition of both quantities being given in the following. 



The hard part of the elementary interaction, CTjj^ddor' i^ given as 
^laddcr(s) = J2 dw+dw-dQl 



(2) 






where da3Qj,j^/(iQ^represents the hardest scattering in the middle of the ladder (indicated symbolically by 
the somewhat thicker ladder rung in fig. ll), and where £^qcd represents the evolution of parton cascade 
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Figure 2: The calculation of -Eqcd- 
from scale Q\ to Q\, using the DGLAP approximation | alt82 , |ell96 |, given as (see fig. ||) 
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where -Eqcd represents an ordered ladder with at most n ladder rungs. This is calculated iteratively based 
on 



^tcl^Ql Ql ^) = -5(1 - x) 5., ^\Ql, Ql) 
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where the indices i, j, k represent parton fiavors. P^ are the Altarelli-Parisi splitting functions and A^ 
is the so-called Sudakov form factor. The soft part of the elementary interaction, Egoit, is the usual soft 
Pomeron expression. 

In addition to the semihard contribution Uscmi, one has to consider the expression representing the 
purely soft contribution: 

-Psoft {s,x~^,x~ ,b) = £;soft ( ^ °_ ,b] , (5) 

with the scale parameter sg = 1 GeV. The complete contribution, representing an elementary inelastic 
interaction in an energy range of say 10 - 10'' GeV, is therefore given as 



D = D 



soft 
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(6) 



We would like to stress, that the "soft end" of the semihard Pomeron has exactly the same structure as 
the soft contribution itself, no new parameters enter. 

There is still something missing: the outer legs of the elementary diagram are not the nucleons, but 
nucleon "constituents", to be more precise quark-antiquark pairs. We call these constituents also "par- 
ticipants" to indicate that they are actively participating in the interaction, in contrast to the remnants, 
which represent the non-participating part of the nucleons. So for each incoming leg, we have an additional 
factor Fpart, which we assume to be of the form impart (a;) = x^"?"". Similarly, for each remnant, we add a 
factor F^cmn{x) = a;^"""°, where the arguments of -Flomnla;) are the momentum fractions of the remnants. 
We define 

where we introduced a factor 2 for later convenience. 

3 Multiple Scattering 

We assume that the dominant diagrams for nucleus- nucleus scattering are those which consist of elementary 
diagrams as discussed in the previous sections, see fig. ||. One easily writes down the corresponding 



^inel -^ 



all cuts 
all diagrams 




Figure 3: The inelastic cross section for nucleus-nucleus scattering. 



formula: each cut Pomeron contributes a factor 2CG, each uncut one a factor {—2CG){the semihard 
Pomeron amplitude is assumed to be imaginary), and each remnant F+ or F^ . So we get 
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where / dTAB represents the integration over transverse coordinat es of pr ojectile and target nucleons with 



the appropriate weight given by the so-called thickness functions |wer93|. The factor Dab is given as 



"---U{7c^''-7c'->}U{ic^^'-ic''^}' 



with 



1 if nucleon is passive 
if nucleon is active 



where an active nucleon participates in at least one elementary interaction, whereas a passive one does 
not. The functions n{k) and T(fc) refer to the projectile and the target nucleons participating in the k^^ 
interaction. We fully account for energy-momentum conservation, which we consider extremely important 
due to the nonplanarity of the diagrams, implying the interactions to occur in parallel. 

The expansion of cTinoi in terms of cut diagrams as given in eq. |7| represents a sum of a large num- 
ber of positive and negative terms, including all kinds of interferences, which excludes any probabilistic 
interpretation. Our strategy consist s there fore of performing partial summations such that the remaining 



terms allow such an interpretation [|agk73|l . So we classify the diagrams according to the cut elementary 
diagrams (real emissions), and then sum over all diagrams of a given class, which amounts to summing 
over uncut elementary diagrams (sum over virtual emissions): 

^incl = J2 JdXdXV 

cut diagrams X> 

= E ^^1 E jdxvy 

cut ladders l^ uncut ladders J 

where V is the mathematical expression corresponding to a diagram as shown in fig. y, appearing in 
formula |7|, and where X and X represents all the light cone momenta of the cut and uncut elementary 
diagrams. The term in brackets {. . .} may be finally interpreted as probability for the corresponding 
"ladder configuration". 

Let us write the formulas explicitely. We have 
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The variables appearing in eq. (^) may be represented by two multi variables: the interaction-type variable 
M which specifies for each of the AB nucleon pairs the type of the interaction (how many cut Pomerons 
of which type occur), and the momentum variable X, already mentoned earlier, which specifies for each 
elementary interaction the momentum fractions. Eq. (^) may thus be written as 

fTiuci =Y, fdxn{M,x). (10) 



where D,{M, X) is the integrand of eq. (||). Both variables, K = {M, X} represent a "ladder configuration" 
and V,{M,X) considered to be the corresponding probability density. 
There are two fundamental problems to be solved: 

• the sum over virtual emmisions has to be performed 

• tools have to be delevopped to deal with the multidimensional probability distribution n{M,X). 

Both are difficult tasks. There is no way to do the summation numericly in case of two heavy nuclei 
A and B: we have AB summation indices, so if we assume 10 terms per index to reach convergence, 
we have to sum over 10"^^ terms! It is also out of question to use Monte Carlo methods, since positive 
and negative terms occur. However, we are able to provide a solution, as discussed later. Concerning 
the multidimensional probability distribution n{M,X), we are going to develop methods well known in 
statistical physics (Markov chain techniques), which we also are going to discuss in detail later. So finally, 
we are able to calculate the probability distribution Q{M,X), and are able to generate (in a Monte Carlo 
fashion) "ladder configurations" (M , K) according to this probability distribution. The next task amounts 
to generate explicitely partons, again based on our master formula eq. g. This will be discussed in next 
section. 

4 Parton Configurations 

In this section, we consider the generation of parton configurations in nucleus-nucleus (including proton- 
proton) scattering for a given ladder configuration, which means, the number of elementary interactions per 
nucleon-nucleon pair is known, as well as the light cone momentum fractions x~^ and x~ of each elementary 
interaction. A parton configuration is specified by the number of partons, their types and momenta. We 
showed earlier that the inelastic cross section may be written as 

Keic 

where the symbol S means S / and where K — {M, X} represents a ladder configuration. The function 
D,{K) is known (see eq. (^) and is interpreted as probability distribution for a ladder configuration K. 
For each individual ladder a term 2CG appears in the formula for ^{K), where 2G itself can be expressed 
in terms of parton configurations, which provides probability distributions for parton configurations, and 
which provides the basis for generating partons. We want to stress that the parton generation is also 
based on the master formula eq. (||), no new elements enter. In the following, we want to discuss in detail 
the generation of parton configurations for an elementary interaction with given light cone momentum 
fractions x~^ and x~and given impact parameter difference b between the corresponding pair of interacting 
nucleons. 

First, we have to specify the type of elementary interaction (soft or semihard). The corresponding 
probabilities are 

Gscmi(s, a;+, x", 6) / G{s, a;+, x" , b) (12) 

and 

Gsottis,x^,x^,b) / G{s,x~^,x~,b) (13) 

respectively. 

Let us now consider a semihard contribution. We obtain the desired probability distributions from the 
explicit expressions for 20801111- For given x'^, x~ , we have 

2G,,nAs,x+,x-,b) (X Jdx+dx^ J /d'fo'E Kofti^,b')El,,{^,b-b')ali^^Jx+x^s) \ , (14) 



with 



(s) = Yj f dw+dw-dQ^ (15) 
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representing the perturbative parton-parton cross section, where both initial partons are taken at the 
virtuality Qo- The integrand {. . .} of eq. ^ serves as probability distribution to generate x^ and x^. 

Knowing the momentum fractions x^ and x^oi the "first partons" of the parton ladder, we can construct 
the complete ladder. To do so, we generalize the definition of the parton-parton cross section criaddcr to 
arbitrary virtualities of the initial partons, we define 

^iaddcr(Q^Q2,s) = ^ /'du;+du;-dg2 (16) 

kl •' 

rjfjkl 

and 

^trA^QlQls) = 5]/'d«;-dQ2 (17) 

k 

representing ladders with ordering of virtualities on both sides (crjg^) or on one side only (o'qj.j). We 
calculate and tabulate cr^g^ and cr^j.^ initially so that we can use them via interpolation to generate 
partons. The generation of partons is done in an iterative fashion based on the following equations: 
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5 Outlook 



So far we presented a consistent and very transparent new approch to calculate parton production in 
nucleus-nucleus (including nucleon-nucleon) scattering. But, unfortunately, the real world consists of 
hadrons, so we still have to deal with the problem of hadronization. This is not so clear. We provide a 
"minimal model" where we simply translate the partons from each individual elementary interaction into 
the language of relativistic strings, the latter ones being decayed using the machinery of relativistic string 
decay. An alternative would be to take our partons as initial condition for a transport treatment of a 
partonic system. We do not want to explore these options any further in this paper. 
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